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Introduction 

The study of e^e^ annihilation has always been one of the most useful means of exploring 
QCD and determining a value for as(M^), the strong coupling constant. Multi-jet rate^ in 
particular enable us to examine the perturbative nature of QCD with long distance effects kept 
comparatively low. It is also important to make sure that the shape variable being investigated 
is both collinear and infra-red safe. Jet rates can satisfy all of these criteria as long as care 
is taken in the choice of the jet clustering algorithm. In this paper we calculate the leading 
and next-to-leading logarithmic contribution to the four-jet rate for e^e~ annihilation using the 
Durham algorithm. (For an explanation of various algorithms and the reasoning behind choosing 
the Durham one see refs. ||l|, ^, ^). We then obtain an expression for the jet rate in terms of a 
dimensionless jet resolution parameter, ycut, which can be considered as a measure of how well 
we are able to resolve two approximately collinear partons. According to the Durham algorithm 
we define ycut = Qq/Q^, where Q ~ ^/s is the scale of the jet-production process and hence the 
cut-off energy scale Qq can be considered to be the energy threshold below which the process 
starts to become non-perturbative. 

In the region of small ycuti"^ 1) the emitted gluons are predominantly soft and collinear 
resulting in the logarithmic enhancement of higher orders 0,0. It is therefore necessary to 
resum them to all orders in Ug to obtain a reliable prediction for the four-jet rate. 

Leading Logarithms and Exponentiation 

First it is important to stress what we are actually calculating in the resummation procedure. 
Using the coherent branching formalism |5|, |6|, , we are able to resum exactly all contributions 
to the shape variable at leading logarithms, LL (a"L^") and next-to-leading logarithms, NLL 
(a^L^"~^) in the perturbative expansion where L = — In(ycnt)- This means that all terms 
sub-leading are not completely reproduced and therefore they are dropped in our calculation. 

The idea behind exponentiation is to increase the domain of applicability of the shape variable 
such that it extends into the region of agL < 1. The result of this procedure is to obtain a closed 
function of the form T {Lgi{asL) + g2{asL)), where gi^OgL) resums all leading-logarithmic con- 
tributions and g2{oisL) resums the next-to-leading ones such that when expanded the whole 
perturbation series is reproduced down to terms of the form a'^L'^, where n < m < 2n. For 
the jet fractions being studied, a simple exponentiation does not arise. It therefore only makes 
sense to calculate the LL and NLL contributions of the perturbation series. 

Calculation 

The primary aim of this paper is to find an analytic expression for the 4-jet rate, Ri{ycut)- The 
most simple way to solve this is to work in terms of a generating function defined by 

oo 

</'^(Q,Qo;^) = E«"^n(Q,Qo) (1) 

n=0 



■I^The n-jet rate is defined as Rn{Q) = (T„-j^t/(^hadrons 



where Rn{Q,Qo) is the probability of finding n-partons of a particular type in the final state 
of a process, p, and u is a jet label to distinguish each of the probabilities. In this case we are 
dealing with e~^e~ — >hadrons, therefore ^^"^ ^ ' = [4>q]'^, where (pq is the generating function for 
a single quark to branch, 

Q 1 

— / dzosizq)^ (---] [(pg{zq,Qo;u) -1]. (2) 

To obtain the n-jet rate, i?„, is simply a matter of differentiating the generating function n times 
at n = 0. The n-jet rate is then 



1 / 5 \" 
RniVcut) = :^ [q~] [(t>qiQ^Qo'^u)Y 



(3) 

u=0 



We find from the application of the coherent branching formalism, the generating function obeys 
the following implicit coupled equations [Q]: 

(pqiQ,Qo;u) =u'^expi2 dqTq{Q,q)[(j)g{q,Qo;u) - 1]\ (4) 



and 

4>g{Q,Qo;u) =u expi / dq{Tg{Q,q)[(l)g{q,Qo;u) -1]-Tf{q)} 
\jQo ' 

x^'^ + uj dqTf{q)exp(^l' dq'{[2Tq{q,q') -Tg{q,q')][cl)g{q',Qo;u) - I] +Tf{q')}^y (5) 
Where the emission probabilities are defined as 

vr q \ q 12/ 
Tfy) = ^^. (8) 

The 2-jet limit is important as the jet rate becomes semi-inclusive and exponentiation holds 
exactly. This gives 
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R2{ycut) = exp\^^{3-L)-Po \^ |, (9) 



where L = ln(l/ycut), a = Ois{Q)/TT and we have used Pq = {UCa — 2Nf)/3. The 3-jet case was 
evaluated in [^ and we proceed in a similar way. 



Firstly we find (3) gives, in the n=4 case 

RiiVcut) = '^R2{ycut){j dqTq{Q,q)Ag{q) I dq'Tg{q,q')Ag{q) 

rQ 



+2R2{ycut) f dqrg{Q,q)Ag{q) T dq' rf{q')Af{q') 



+R-i{ycut)\j dqTg{Q,q)Ag{q) 
where we have introduced the Sudakov form factors 
exp I 



AgiQ) = exp - / dqrgiQ,q) 
\ JQo 



AgiQ) = exp(-J^^dq[Tg{Q,q)+Tf{q)]\, 

AfiQ) = exp(-j\[2Tg{Q,q)-rg{Q,q)-Tf{q)]\. 
We need only work with the one-loop definition of the strong coupling constant, 



as{Q) 



1 + Ai^ln(2) 



(10) 

(11) 
(12) 

(13) 

(14) 



Riiy, 



cut) 



Ri 



+ ^07^ 

/3o=0 Opo 



(15) 



as higher order corrections will be sub-leading. Even at this order we are still faced with an 
extremely complicated set of nested integrals. Therefore, as in [g] we proceed by expressing 
Riiycut) as 

SRa 

/3o=0 

This is permissable for any jet multiplicity evaluated at next-to-leading logarithmic order because 
in general we will have 

Rn ^ '^Yi<XLi -\- 'WiiClL -|- • • • 



n r2n— 1 



+ '©n2iiQ' -ty +'©n2n-l'3^ ^ + 



(16) 



where the coefficients ^p, are either j3q independent ("^^2^) or contain a single /3o ('^p2p-i)- All 
other /3o dependence is contained in the strong coupling constant. We note that 



d[as{Q)Y 



dpo 



-m- 



[asiQ)] 



m+l 



/3o=0 



27r 



.ln(^)~a-+iL. 



(17) 



It is now apparent that beyond the first derivative there will only be terms of the form a"-L^" ^ 
which in the NLL approximation can be dropped. The assumption of (15) then seems to be a 



valid one. In fact this expansion greatly simplifies the calculation by enabling us to work with 
terms evaluated with /3o equal to zero. In doing so the coupling Og can really be treated as a 
constant and hence no longer depends on the integration variable. Proceeding in this way, we 
then calculate the four jet rate to be 

R^ = g^(e-2(^+^)((e^-l)'(2 + 3C^aL) 

-x/^e^(e^ - l)(-3 + 12F + 2L)V^erf(\/I))) 
+§(^e~2^( - 24e-2^(e^ - 1)(2 + SaCpL) 

-4y^e~^(2 + 3e^(-3 + 12F + 2L))V^erf(\/I) 
-(12 + aL{llCA - 6Cf{-9 + 12F + 2L)))7r erf(\/I)^ 
+2y^(-7 + 72F + 12L)V2^ evi{V2A))^ 'T^ 

+f (i^l^(e"'^^^^H - 2v/^e^(3 + 2A{-3 + 4F))V^erf(^/I)) 

(1 + e^( - 1 + VlV^erf(VI)))) 
+g:_^(e'2(^+^)2Cf aL(l + 2e^ + 4(e^ - 1)F) 
(1 + e^( - 1 + ^/IV^erf(^/I)))) 

+A{-4 - 2e^(3 - 8F)))0Ferf(VI) 
-e2^(9 - 4^(3 - 8F))\/2^erf(\/2A))) 
+3;^CpaLe"2(^+^)(l - e^(l - 3e^) - 8F(1 - e^) 

+2Fe2^7rerf2(yi)) + (^j, 

erf(x) is the error function defined to be -y= J^ e~y dy and erfi(x) = erf(ix)/i. We have also 

defined A = Caul's/A, F = CpaL'^/A and Q{x,z) = x J^ e^^ erf(y)dy. Attempts were made to 
solve Q exactly but no closed form was found. It appears that the integral is just a generalisation 
of the error function and hence cannot be solved, except in certain cases. A reference containing 
various properties of this function is given [0]. 



Properties of the Four-Jet Rate 

With the complete result calculated we are able to reproduce the exact LL and NLL coefficients 
of Us at any order. The first three orders in Os/vr are given. 

i.e. RAiVcut) = a^B^L^ + B^L^ + 0{L^)) + 03(^6^^ + C^L^ + 0{L^)) + a^DgL^ + DrL'^ + 
0(L6)) + .... 



B4 = gCp + j^CfCa- 

-B3 = ^Cp — j^CfCa + -^CpNf. 

r~i —1/^ 3 1 /^ 2/^^ 7 /'^ /'^ 2 

(_/6 — le'L'F — As^F ^A — 2880^^*-"^ • 

C5 = tbCf' + Ts^F Ca + ^CfGi - ^Cf Nf - j^CACpNf. 

-^8 = q^Cf + TasCp Gl + 5J2CF Gl + 5T2o'^^Gl • 

n _ -Sr* 4 17^ 3ry 1439 r< 2^2 2371 /"y /-< 3 , 323 /^ 3 Ar i 31 r» 2^^ AT 

^7 - 16<^F - ejCvF 04 - 17280*-^^ '-^^ ~ 241920^^^^ + T0080^^ ^^/ + 3024^^ ^^^^/ 



This is in agreement with ||l| which gives the -64^3 coefficients. The 6*6,5 coefficients were in 
addition calculated by expanding out the integral equation (10) as a function of Og. Another 
test was to calculate R4 in the large Nc limit {Nc is the number of colours) to the order of 
leading logarithms. This greatly simplifies the equations as Gi — > Nc, Cf — > Nc/2 and Nj can 
be disregarded. Eq. (4) now collapses down to 



2 / /"^ 
4>{Q,Qo',u) = u exp 2 / dgr„((5 

V JQo 



-(j){q,Qo;u) - 1 
u 



(18) 



Also noting that at leading logarithmic order R^ will be independent of Pq we can safely set it 
to zero. We then get 

i?f- = ie-3^ (e - 8e^ + 4VIe^(l - 2e^)V^eriVA - (1 - 2A) e^^ erf ^^/I 

+26^^(1 + 2\/I\/2^erf\/2]4)) . (19) 



This is in agreement with the full NLL result in the appropriate limit. 

We also note that in the psuedo-abelian limit of simply Ca and Nj going to zero, exact 
exponentiation holds. We also find that this gives a reasonably good approximation to the full 
non-abelian case within about 15-20%. 



Conclusion 

To conclude, we have found an analytic expression for R4,{ycut) which exactly resums the leading 
and next-to-leading logarithmic contributions to all orders in Ug- Previous work on including 



the leading and next-to-leading logarithms to all orders [10| was performed by solving (10) 
numerically. We propose that it is beneficial to have an explicit result since it is then possible to 
extract the purely LL and NLL contributions and drop the incomplete sub-leading terms. Most 
importantly we will be able to utilise the result to address the so called 'renormalisation scale 
ambiguity' in combining the resummed result with the fixed order one. We will be implementing 
a novel technique and hope to report on this in the future. 
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